function [v] = get_beads_volume(LF, x, y, z, sizexyz)
% x , y, z : beads loc
% sizexyz : beas size

[sizey, sizex, sizez] = size(LF);
half = floor(sizexyz/2);

if isempty(find([x, y, z] + half > [sizex, sizey, sizez])) && isempty(find([x, y, z] - half < [1,1,1]))
    v = LF(y-half(2):y+half(2), x-half(1):x+half(1), z-half(3):z+half(3));
    v = double(v);
    %     v = v/max(v(:));
else
    v = 0;
end

end

